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Abstract 

The data taken by Delphi at centre-of-mass energies between 189 and 209 GeV 
are used to place limits on the CP-conserving trilinear gauge boson couplings 
Agf , A 7 and Ak 7 associated to W + W~ and single W production at Lep2. 
Using data from the jjiv, jjjj, jjX and £X final states, where j, £ and X 
represent a jet, a lepton and missing four-momentum, respectively, the following 
limits are set on the couplings when one parameter is allowed to vary and the 
others are set to their Standard Model values of zero: 

Agf = -0.025™ , 
A 7 = 0.002+H! 
and Ak 7 = 0.024±g;gS . 

Results are also presented when two or three parameters are allowed to vary. 
All observations are consistent with the predictions of the Standard Model and 
supersede the previous results on these gauge coupling parameters published 
by Delphi. 



This paper is dedicated to the memory of Professor Paul Booth who was Delphi 
Deputy Spokesperson from 1989 to 1994- He played a key role in the final 
installation and commissioning of the Delphi detector as well as leading the 
Liverpool group for many years thereafter. 
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1 Introduction 

The reactions e + e~ — > W + W~ and e + e~ — > Weu can be used to test the non-Abelian 
nature of the Standard Model (SM) by studying the trilinear couplings of the electroweak 
bosons [1]. In this paper, data from the final states jjiu, jjjj, jjX and £X (where j 
represents a quark jet, £ an identified lepton and X missing four-momentum) taken by 
Delphi at centre-of-mass energies from 189 to 209 GeV are used to determine the values 
of the coupling parameters which describe the trilinear WW Z and WW^y interactions. 

The WWV vertex (V = Z or 7) can be described by an effective Lagrangian with 
14 parameters [1,2]. The set of parameters is reduced to five by assuming electromag- 
netic gauge invariance and by restricting the contributions in the effective Lagrangian to 
operators which are C, P-conserving. A further reduction is then achieved by extract- 
ing from the CP-conserving Lagrangian those terms which satisfy SU(2) ® U(l) gauge 
invariance, are not constrained by existing low-energy data, and are of lowest dimen- 
sion (< 6). This leads to a set of three independent parameters, which are studied by 
Delphi in the present paper: Agf , the difference between the overall WW Z coupling 
and its SM value, A/t 7 , the deviation of the dipole coupling « 7 from its SM value, and the 
quadrupole coupling, A 7 . The imposition of gauge invariance implies relations between 
the dipole couplings « 7 and Kz and between the quadrupole couplings A 7 and Xz, namely: 
Ak z = Agf — rc" 2 g^ A/t 7 and \ z = A 7 , where 9 W is the electroweak mixing angle. The 
terms in the effective Lagrangian which conserve CP, as well as C and P separately, 
correspond to the lowest order terms in a multipole expansion of W-^ interactions: 



Qw = egj , (1) 

g 

hw = t, — (gl + ^ + Xj) (2) 

2m w 

g 

and qw = — — A 7 ) , (3) 

where Qw, l^w, and qw are respectively the charge, the magnetic dipole moment, and 
the electric quadrupole moment of the W + . It may be noted that electromagnetic gauge 
invariance, invoked above, implies the value g\ = 1 in these relations. 

The diagrams which contribute to W + W~ production are shown in figures 1(a) 
and (b). The WWV vertex only occurs via the s-channel diagram shown in figure 1(a) 
and not in the t-channel diagram, shown in figure 1(b), which leads to the same final 
states. This reaction is studied in this paper in the final states where one W boson de- 
cays to hadrons and the other decays into leptons, jj£v, and when both W bosons decay 
into hadrons, jjjj. The WW^y vertex alone is also accessible at Lep2 through single 
W production and is shown in figure 1(c). This process contributes significantly in the 
kinematic region where the final state electron is emitted at a small angle and is studied 
here in two final state topologies: £X, where the W boson decays into a lepton and a 
neutrino, and jjX, where the W decays into a pair of quarks. 

Delphi has previously published results on charged trilinear gauge coupling parame- 
ters using data from WW and Wev production at energies up to 189 GeV [3-5], and a spin 
density matrix analysis of Delphi data from the jjev and jj\xv final states at energies 
up to 209 GeV has been used to determine both CP-conserving and CP-violating cou- 
plings [6]. The results presented here supersede all those on CP-conserving couplings in 
these publications. Results at energies up to 209 GeV from the other Lep collaborations 
can be found in references [7-10]. 

The Delphi detector is described in section 2. The data and simulation samples 
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Figure 1: Diagrams contributing to W + W~ and W^ez/ production at Lep2. (a) and (b) 
are the diagrams which describe W + W~ production and (c) describes Weu production. 
The trilinear gauge boson vertices are denoted by shaded circles. 

describes the analysis techniques used in the extraction of the couplings from the data. 
The different sources of systematic uncertainties are discussed in section 6 and the results 
from fits to the data are given in section 7. The conclusions are presented in section 8. 



The Delphi detector and its performance are described in detail in [11,12]. For Lep2 
operation a number of changes were made to the sub-detectors, the trigger [13], the run 
control system and the track reconstruction algorithms to improve the performance. The 
angular coverage of the Vertex Detector was extended [14] to cover polar angles 1 in the 
range 11° < 9 < 169° with the inclusion of the Very Forward Tracker. Together with 
improved tracking algorithms, alignment and calibration procedures, this resulted in an 
increased track reconstruction efficiency in the forward region of Delphi. 

During the final year of operation, one sector of the twelve that constituted the central 
tracking device (TPC) ceased to function. This affected around a quarter of the data 
collected in 2000. The tracking algorithms were modified in this sector so as to reconstruct 
tracks from the signals in the other tracking detectors. 

3 Data samples 

A total integrated luminosity of around 600 pb _1 was collected by Delphi between 
1998 and 2000. Table 1 shows the integrated luminosity available at each energy and 
the luminosity-weighted centre-of-mass energies. The luminosity was determined from 
Bhabha scattering measurements [15]. 

All four-fermion final states were generated with the four-fermion generator 
WPHACT [16,17], set up as described in [18]. The most recent radiative corrections 
to the W pair production cross-section, calculated in the so-called Double Pole Approxi- 
mation (DPA), were included via an interface to YFSWW [19]. 

The background from two-fermion production was simulated using KK2f [20] and 
KoralZ [21]. Additional background contributions from two-photon production were gen- 
erated using BDK [22] and BDKRC [23]. All of the generators were interfaced to the 
PYTHIA [24,25] hadronisation model tuned to the Delphi data collected at the Z res- 
onance [26]. 

1 The Delphi coordinate system has z axis in the direction of the incoming e~ beam. The polar angle 9 is defined with 
respect to this direction, and the r<j> plane is perpendicular to the z axis. 



2 The DELPHI detector 
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Luminosity-weighted y/s (GeV) Hadronic L (pb x ) Leptonic L (pb x ) 



188.63 


154.4 


153.8 


191.58 


25.2 


24.5 


195.51 


76.1 


72.0 


199.51 


82.8 


81.8 


201.64 


40.3 


39.7 


204.81 


82.6 


74.9 


206.55 


135.8 


123.7 



Table 1: The centre-of-mass energies weighted by the integrated luminosity (L) for each 
of the Lep2 data taking periods. The hadronic luminosity was used in the fully hadronic 
selection and the leptonic luminosity was used in the other channels. The different lu- 
minosities are due to tighter requirements being made on the detectors used for lepton 
identification in the semi-leptonic channels. 



The large simulated samples (about 1M charged current four-fermion events, 500K 
neutral current four-fermion events and 1M two-fermion events at each energy) were 
interfaced to the full Delphi simulation program DELSIM [11,12] and passed through 
the same reconstruction chain as the experimental data. In order to allow analysis of 
data taken during the period when one part of the TPC was inoperative (as described in 
section 2 above), additional samples were generated with the detector simulation modified 
to model this situation. 



4 Event selection 

In this section, the selection of events in the various final state topologies used in the 
determination of the coupling parameters is described. 

Events selected for analysis from the WW final state came from the semi-leptonic 
channel, jjlv, and from the fully hadronic channel, jjjj. The semi-leptonic final state was 
divided into three further channels, jjeu, jj/jtv and jjTu; in the case of jjrv production, 
only events with the tau decaying into a single charged track were considered. Events in 
the semi-leptonic final state are therefore characterised by two or more hadronic jets, an 
isolated lepton - this comes directly from the decay of the W or from the cascade decay 
of the tau lepton - or a low multiplicity jet due to a hadronic tau decay, and missing 
momentum from the neutrino(s). The main backgrounds come from qq(^y) production 
and from four-fermion final states of two quarks and two leptons of the same flavour. 
The analysis of the fully hadronic final state from WW production involved a search for 
four-jet events in which the di-jet invariant masses of one of the pairings into two di-jets 
were compatible with the W mass. Here, also, 55(7) production represents a major source 
of background, with some contamination also from ZZ decays into qqqq and qqr + r~ . 

Events from single W production, Weu, were selected in the kinematic region where 
the final state electron is very close to the beam direction and remains undetected. Final 
states with hadronic W decays and with leptonic decays into electron or muon and 
a neutrino were considered, so that the topologies analysed were jjX and £X, with 
I = e, n and X representing missing momentum. The main background contributing to 
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the jjX topology came from qqrv production. In the eX channel the major source of 
background was from e + e~7(7) production with one electron (or positron) and the final 
state photon(s) unobserved, while in the \iX topology the main backgrounds were from 
ee/ji/j, production, mainly via two-photon processes, and from fifi'-f production. Some of 
these background processes (such as qqrv production) themselves contain triple gauge 
boson vertices in their production mechanisms, and thus contribute to the precision of 
the results. 

Full details of the event reconstruction procedure adopted, and of the selection of 
events in the channels considered here from WW production can be found in [27], 
Delphi's report on the measurement of the WW production cross-section, while the 
selection procedure for events in the Wev final state is very similar to that used in our 
previous publication of charged trilinear gauge boson couplings at 189 GeV [5]. In the 
following sections, a summary of these procedures is given. 

The total numbers of events selected at each centre-of-mass energy are given in table 2. 
The table also gives examples (at 200 GeV) of the event selection efficiencies and estimated 
background cross-sections; the errors on these cross-sections are treated as a systematic 
uncertainty and are discussed in section 6. 

4.1 Particle selection 

Reconstructed charged particles were required to have momentum greater than 
0.1 GeV/c and less than 1.5 times the beam momentum, a relative momentum error 
less than 1, an impact parameter in rip less than 4 cm, and a z impact parameter less 
than 4 cm/sin Neutral clusters were required to have energy exceeding 300 MeV in 
the barrel electromagnetic calorimeter (HPC) and exceeding 400 MeV and 300 MeV in 
the two forward electromagnetic calorimeters (FEMC and STIC, respectively). Elec- 
tron identification was based on the association of energy deposits in the electromagnetic 
calorimeters with momentum measurements in the tracking chambers (the Inner Detector 
and the TPC) and, in the case of lower energy candidates, with energy loss measurements 
in the TPC. Muon candidates were identified by extrapolating tracks through the entire 
detector and associating them with energy deposits recorded in the hadron calorimeter 
(HCAL) and hits recorded in the muon chambers. 

4.2 Selection of events in the jjiu final state 

The selection of events in the semi-leptonic final state involved cut-based selections, 
followed by the application of an Iterative Discriminant Analysis (IDA) [28,29]. 

An initial hadronic pre-selection was applied where at least 5 charged particles were 
required, the energy of the charged particles had to be at least 10% of the centre-of-mass 



energy, and the following condition was imposed: ^jEMF'j + EMF£ < 0.9xE beam , where 

EMFf y b are the total energy deposited in the electromagnetic calorimeters in the forward 
and backward directions, defined as two cones around the beam axes of half-angle 20° . 

At this point a search was made for leptons, allowing each event to have up to three 
lepton candidates, one of each flavour. Of the electrons found in the particle selection 
procedure, the one with the highest value of E x 9 iso was chosen as the electron candidate. 
Here E is the measured electron energy and 6 iso is the isolation angle of the electron track, 
defined as the angle made to the closest charged particle with momentum greater than 
1 GeV/c. The candidate was then required to have energy greater than 15 GeV. Similarly, 
an identified muon track with momentum p was selected if it had the highest value of 
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p x 6 iso and if its momentum exceeded 15 GeV/c. The event was then clustered into jets 
using LUCLUS [24,25] with dj oin = 6.5 GeV/c. Particles were removed from the jets if 
they were at an angle greater than 20° to the highest energy particle and the remaining jet 
with the lowest momentum-weighted spread 2 was considered as a tau candidate. Particles 
were removed from this jet if they were at angle greater than 8° from the jet axis and the 
remaining jet was required to contain at least one charged particle. 

For each lepton candidate, the remaining particles in the event were clustered into two 
jets using the DURHAM algorithm [30]. Each of these jets was required to contain at 
least three particles, of which at least one had to be charged. A further pre-selection was 
made before applying the full selection using the IDA: for jjev and jjfiv candidates the 
transverse energy was required to be greater than 45 GeV; the missing momentum had 
to exceed 10 GeV/c; the visible energy divided by the centre-of-mass energy at which the 
IDA was trained (defined below), E vis / E train , was required to be between 40% and 110%; 
and the fitted W mass from a constrained kinematic fit (imposing four-momentum con- 
servation and equal mass for the two W bosons in the event) had to be greater than 
50 GeV/c 2 . For jjrv candidates, the transverse energy was required to be greater than 
40 GeV, the missing momentum between 10 and 80 GeV/c, the ratio E vis /E train between 
35% and 100%, and the fitted W mass greater than 50 GeV/c 2 . 

After the pre-selection cuts an extended IDA analysis was used which treated correctly 
the case where the signal and background had different shapes. The input observables 
were transformed to make their distributions Gaussian. The IDA was trained on 50k four- 
fermion events for charged and neutral processes and 100k qq^) events at three centre-of- 
mass energies: 189, 200 and 206 GeV. The following variables were used in the selection 
of all channels: the total multiplicity, the visible energy, the lepton isolation angle, the 
ratio between the reconstructed effective centre-of-mass energy, y/s', [31] and the centre- 
of-mass energy, y/s, the magnitude and the polar angle of the missing momentum, and the 
fitted W mass. The lepton energy was used in the selection of jjev and jj\iv candidates. 
The angle between the lepton and the missing momentum was used in the jj/iv and jjrv 
selections. For the jjev selection, the transverse energy was also used. In addition, the 
aplanarity 3 , the charged multiplicity of the tau jet and its momentum-weighted spread 
were used in the jjrv selection. The cut on the output of the IDA was chosen such that 
it maximised the value of the efficiency times the purity for each channel. 

In the application of the IDA, events with more than one lepton candidate, one of 
which was a muon, were first passed through the jj/iv selection procedure; those not 
selected, but containing an electron candidate, were then passed to the jjev selection, 
and if the event failed both the muon and electron selection procedures and included 
a tau candidate, it was passed to the jjrv selection. A final cut, requiring the charged 
multiplicity of the tau jet to be 1, ensured that the charge of the W boson which produced 
it was well determined. 

At centre-of-mass energy of 200 GeV, the efficiencies of the jjev, jjfiv and jjrv 
selections were found to be 71.0%, 88.2% and 54.6%, respectively (see table 2). The 
selection efficiencies differed by no more than 2% over the energy range considered. The 
respective background cross-sections for the three channels at 200 GeV were evaluated to 
be 0.232 pb, 0.075 pb and 0.344 pb, with the main contributions coming from qqi^j) and 
from neutral current four-fermion final states. Combining the data at all centre-of-mass 

2 defined as ^A^' i P | ' where 0; is the angle made by the momentum Vi of the ith particle in the jet with the total jet 
momentum. 

3 defincd as |A3 where A3 is the smallest eigenvalue of the sphericity tensor S a P = ^ A^-A ■ The p; are the three 
momenta of the particles in the event and a, f3 = 1, 2, 3 correspond to the x, y, z momentum components. 
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energies, totals of 1101, 1246 and 886 events were selected in the three leptonic channels, 
respectively. 

4.3 Selection of events in the fully hadronic final state 

In the selection of fully hadronic final states, the charged and neutral particles in 
each event were forced into a four-jet configuration with the DURHAM algorithm. A 
pre-selection was performed where the reconstructed effective centre-of-mass energy, y/s', 
was required to be greater than 65% of the nominal centre-of-mass energy, the total and 
transverse energy for charged particles were each required to be greater than 20% of the 
nominal centre-of-mass energy, the total multiplicity for each jet had to exceed 3, the 
condition y cut > 0.0006 was imposed for the migration of 4 jets to 3 jets when clustering 
with the DURHAM algorithm, and a four-constraint kinematic fit of the measured jet 
energies and directions, imposing four-momentum conservation, was required to converge. 

A feed-forward neural network, based on the JETNET package [32] was then used 
to improve the rejection of two- and four-fermion backgrounds. The network uses the 
standard back-propagation algorithm and consists of three layers with 13 input nodes, 
7 hidden nodes and one output node. The choice of input variables was optimised [33] to 
give the greatest separation between WW and two-fermion events. The following jet and 
event observables were used as input variables: the difference between the minimum and 
maximum jet energies after the kinematic fit, the minimum angle between the jets after 
the fit, the value of y cut from the DURHAM algorithm for the migration of 4 jets into 
3 jets, the minimum particle multiplicity of any jet, the reconstructed effective centre- 
of-mass energy, the maximum probability amongst each of the 3 possible jet pairings of 
a six-constraint fit (imposing the additional constraints that the invariant mass of each 
jet pair should be equal to the W mass, set equal to 80.40 GeV/c 2 ), the thrust, the 
sphericity, the transverse energy, the sum of the cubes of the magnitudes of the momenta 
of the 7 highest momentum particles, Y^l=i Ip«| 3 > the minimum jet broadening, B min [30], 
and the Fox- Wolfram moments if 3 and HA [34]. 

The neural network was trained on separate samples of 2500 signal and Z/7 — > qq 
events for each centre-of-mass energy. The network output was calculated for other 
independent four-fermion, two-fermion and two-photon processes. 

The efficiency of the fully hadronic selection for a centre-of-mass energy of 200 GeV 
was estimated to be 81.9% (see table 2); the efficiency varied by no more than 4% over 
the energy range considered. The background cross-section at 200 GeV was evaluated to 
be 1.21 pb, with the main contribution coming from qq^). Combining the data at all 
centre-of-mass energies, a total of 4348 events was selected. 

4.4 Selection of events in the jjX final state 

Events were considered as jjX candidates if there were no identified leptons with mo- 
mentum greater than 12 GeV/c, the measured transverse momentum exceeded 20 GeV/c, 
and the invariant mass of detected particles lay between 45 GeV/c 2 and 90 GeV/c 2 . In 
addition, events were rejected if any neutral clusters were found in the electromagnetic 
or hadronic calorimeters with energy exceeding 1 GeV within a cone of half-angle 30° 
around the direction of the missing momentum. Particles were clustered into jets using 
LUCLUS with dj oin = 6.5 GeV/c and events were required to have two or three jets only. 
Surviving events were forced into a two-jet configuration and accepted if the jet polar 
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angles were between 20° and 160° and the acoplanarity angle 4 between the jets was less 
than 160°. 

The efficiency of the selection is quoted with respect to a reduced phase space defined 
by the following generator level cuts: the acoplanarity angle between the quarks was 
required to be less than 170°; the invariant mass of the quark pair had to be greater 
than 40 GeV/c 2 ; the quark directions were required to have polar angles between 20° and 
160°; and the electron polar angle was required to be less than 11° or greater than 169°. 
The efficiency for selecting the Wev final state with W— > qq was found to be between 
43.7% and 48.0%, depending on the centre-of-mass energy, with a luminosity-weighted 
mean value of 45.4%; 215 events in total were selected in the data. For Standard Model 
values of the couplings, a total of 219.8 ± 1.6 events was expected, comprising 79.5 events 
from qqev production with the electron or positron lost in the beam pipe, 13.0 events 
from qqev production with the electron or positron elsewhere in the detector, 18.5 events 
from qqfiu production, 67.0 events from qqrv, 36.0 events from qqvv, and 5.8 events 
from qql'-f) production. The error in the expected total number of events arises from 
the statistical errors in the selection efficiencies estimated for the contributing processes. 
All the processes contributing to the selected sample except qq^) production include 
diagrams with trilinear gauge couplings, and this was taken into account in the subsequent 
analysis. The background cross-section of 0.048 pb shown for the jjX channel in table 2 
represents the contribution at y/s = 200 GeV from qqij) production. 

4.5 Selection of events in the £X final state 

To be considered as an tX candidate, events were required to have only one charged 
particle, clearly identified as an electron or a muon from signals in the electromagnetic 
calorimeters or the muon chambers, respectively, using the same procedures as described 
in the selection of semi-leptonic events (section 4.2). The impact parameter for the 
lepton was required to be less than 0.1 cm in the r<p plane and less than 4 cm in the z 
direction. The lepton candidate was required to have momentum less than 75 GeV/ c, with 
the transverse component of this momentum greater than 20 GeV/c. The total energy 
deposited in the electromagnetic calorimeter, but not associated with the track, was 
required to be less than 5 GeV. The ratio of the energy deposited by electron candidates 
in the electromagnetic calorimeter to that determined from the measured value of the 
momentum was required to exceed 0.7. 

As in the case of the jjX final state described above, the efficiency of the selection 
was calculated in a reduced phase space region defined by cuts made at generator level; 
for the tX final state, these were defined as follows: the lepton energy was required to 
be less than 75 GeV; the transverse momentum of the lepton had to be greater than 
20 GeV/c; and the polar angle of the missing momentum was required to be in the range 
from zero to 11° or between 169° and 180°. Totals of 37 and 39 candidates were selected 
in the [iX and eX channels, respectively, with luminosity-weighted average efficiencies 
for selection of the Wev final states of 49.4% for W-^-fiv and 31.3% for W^tev. For 
Standard Model values of the couplings, 34.0 ± 1.4 and 31.9 ± 1.5 events were expected 
in the two channels, respectively. The predicted fiX sample comprised 17.7 events from 
e\ivv production with the electron or positron lost in the beam pipe, 1.6 events from ervv 
production, also with an invisible electron or positron, 1.9 events from [i[ivv production, 
2.0 events from jirvv production, 4.0 events from /z/zee, and 6.8 events from fifi^). In 
the eX sample, 19.2 events were expected from eevv production with one lost electron or 

4 dcfincd as the angle between the planes containing each jet direction and the beam direction. 
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Energy (GeV) jjeu jjfxu jjru jjjj jjX jiX eX 



1 8Q 




ooo 


ZOO 


1 D/19 




1 1 
1 1 


i n 

1U 


192 


42 


53 


37 


187 


4 


1 


1 


196 


151 


166 


116 


532 


22 


6 


5 


200 


162 


190 


145 


614 


24 


6 


6 


202 


94 


89 


57 


317 


12 


5 


3 


205 


169 


153 


94 


657 


89 




14 












8 


207 


214 


259 


201 


999 








TOTAL 


1101 


1246 


886 


4348 


215 


37 


39 



e(y/s = 200 GeV) (%) 71.0 88.2 54.6 81.9 48.0 50.8 31.9 
(?back(Vs = 200 GeV) (pb) 0.232 0.075 0.344 1.21 0.048 0.016 0.013 

Table 2: The numbers of events selected from the data in each channel for each centre-of- 
mass energy. The selection efficiencies e and background cross-sections aback are shown 
for the centre-of-mass energy of 200 GeV. 

positron, 1.6 events and 0.7 events, respectively, from e/iuP and eruu with the electron or 
positron in the beam pipe, 1.5 events from ervv production with the electron or positron 
elsewhere in the detector, and 8.8 events from Compton and Bhabha scattering with 
only one electron (or positron) detected in the final state. The background cross-sections 
in the /iX and eX final states at 200 GeV quoted in table 2, 0.016 pb and 0.013 pb, 
respectively, represent the contributions from the processes contributing to these final 
states which have no dependence on the trilinear gauge couplings under consideration, 
namely the /x/x(7) contribution to /iX and the Compton and Bhabha contributions to 
eX. All the other contributions to these final states have a dependence on trilinear gauge 
couplings in their production, and this was taken into account in the subsequent analysis. 



5 Determination of the couplings 

The extraction of the couplings from the data exploited the fact that the differential 
cross-section, is quadratic in the set of couplings, a,(= A</f , Ak 7 , A 7 ), and can be 
expressed as 

^ = cd(n) + + E > ( 4 ) 

where Q represents the kinematic phase space variables and i, j are summed over the num- 
ber, N, of parameters being determined. The coefficients c\ and were calculated using 
WPHACT for the final states coming from W + W~ production and using DELTGC [35] 
for single W final states. This allows the fully simulated events to be re-weighted to 
non-SM values of the couplings. 
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5.1 Semi-leptonic final state 

The analysis of the data in the semi-leptonic channel used the method of Optimal 
Observables [36-38], in which an expansion of the form (4) represents the first two terms 
in a Taylor expansion of the differential cross-section for any process in terms of a set of 
N parameters aij. If it is known that the ctj are small, then the N lowest order terms 
in (4) contain most of the information needed for the determination of the parameters. 
In the present case, where the amplitude for the processes we consider is linear in the 
parameters, the Taylor expansion is truncated at the second order, and (4) gives the 
value of the cross-section without approximation. This suggests an analysis in terms of 
the quantities 

u\(Q) = ^B- and = ^B- , (5) 

c (Q) c (fi) 

which are easily derived from the differential cross-section. Such an analysis is described 
in [39], where the probability distribution function, P(Q,a), for observing an event at 
phase space position Q when the parameters have values a (= a± ... ajv) is projected in 
the u\(£l) and ^(Vt) of (5), the Optimal Variables. 

— * 

When fl is known precisely, a fit to the Optimal Variables allows the couplings to be 
determined with a precision equal to that of an unbinned maximum likelihood fit over all 
of the phase space variables. In practice, the measured values of the Optimal Variables 
are defined by the convolution of the differential cross-section with the resolution and 
efficiency functions of the detector. However, it has been confirmed by Monte Carlo 
tests [40] that little loss of precision occurs when this convolution is performed. In the 
case where one parameter, a iy is free to deviate from its Standard Model value, two 
Optimal Variables (u>\ and ou™ ) contain the whole information, but five (or nine) Optimal 
Variables are required when two (or three) parameters are released from their Standard 
Model values. For one-parameter fits, there is an obvious advantage in simplicity in the 
use of the Optimal Variable method over an analysis using the five angular variables (the 
W production angle and the W + and W~ decay angles) known to contain most of the 
information on the coupling parameters in WW production, while in the case of multi- 
parameter fits the number of Optimal Variables is equal to or greater than the number 
of angular variables. We have compared these methods using simulated event samples: 
in all cases - for one-, two- and three-parameter fits - the precision obtained from the 
Optimal Variable analysis was at least as good as that from the angular analysis, allowing 
us to use the same methodology throughout the analysis. 

The distributions of the Optimal Variables used in fits to the parameters Agf {i — 1), 
A 7 (i = 2) and Ak 7 {i = 3) are shown for the real data and for events simulated with 
SM and non-SM values of the couplings in figures 2 to 5 for a centre-of-mass energy of 
200 GeV. 

The values of the coupling parameters were determined by binned extended maximum 
likelihood fits to the relevant Optimal Variables. A clustering technique was used to 
define the binning of the data, full details of which can be found in [41]. The method 
used the data points to divide the phase space into equiprobable, multidimensional bins. 
For each fit, a set of d variables (d — 2, 5 or 9, as described above) was required to describe 
an event completely and for n events the clustering technique divides the d- dimensional 
space into n bins, each centred on one data point. The available simulated events are 
then assigned to the bins by calculating the scalar distance D M of each simulated event 
k to each of the data points /, 
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D kl = (R t - r k ) M(^-rl), (6) 

and assigning the kth simulated event to the bin / for which D^i is a minimum. In (6), R 
and r are the ci-dimensional vectors that describe the real data point and the simulated 
event, respectively, and M is a d x d matrix representing the metric of the space. The 
metric M was defined by the variances and correlations of the real data distributions of 
the Optimal Variables being determined in any particular fit, so as to take into account 
the fact (observed in figures 2 to 5) that the different variables span different numerical 
ranges. 

The technique described above assumes that the phase space variables Q are fully 
determined for each event. In fact, one ambiguity remains for every event, namely that 
it is not known which of the jets from the hadronic W decay comes from the quark, and 
which from the antiquark. In the analysis, each event was therefore entered twice into 
the maximum likelihood function, once with each of these two assignments. 

A second analysis was performed in the jjiu channel as a cross-check. In this analysis 
a binned maximum likelihood fit was made to the differential cross-section of three angles: 
cos6'vy+, the W + production angle, cos 9i, the polar angle of the lepton with respect to 
the incoming e ± of the opposite sign, and cos^v^; the cosine of the angle between the 
hadronic W and the lepton. The same event selection criteria were applied and the same 
re-weighting method was used as in the Optimal Variable analysis. The distributions of 
these angular variables are shown in figures 6 to 8 for the data and for events simulated 
with different values of the couplings at 200 GeV. 



5.2 Fully hadronic final state 

The analysis of events in the fully hadronic state is complicated by the fact that the 
four observed hadronic jets cannot immediately be assigned to a particular W + or W~ 
decay. Two problems arise from this feature, first, that it is not clear which of the three 
possible pairings of the four jets corresponds to a WW pair, and, second, once the pairing 
is decided, which of the di-jet pairs is the W + and which the W~ . 

The first of these problems was approached by forcing the selected events into a four- 
jet configuration and constructing a neural network to determine the combination which 
was most likely to represent a W pair event. A kinematic fit, imposing four- momentum 
conservation and equal mass for the two di-jet pairs, was performed for each of the three 
combinations. The \ 2 of the kinematic fit and the difference between the nominal W mass 
and the di-jet mass from the fit were used as inputs to the neural network to choose the 
most likely combination. The efficiency of this procedure was estimated to be about 79%, 
where the uncertainty in the pairing was estimated by repeating the procedure using 
simulated events generated with the different parton shower and fragmentation models 
implemented in PYTHIA, HERWIG [42] and ARIADNE [43]. 

The second problem - to distinguish which of the di-jet pairs came from the W + and 
W~ - was partly resolved by constructing an effective jet charge Qj et from the charge of 
the particles in the jet, weighted by their momentum: 

Ql " - e,(p. • JW» ' ( ' 

where and Pi are, respectively, the charge and momentum of the particle in the jet, 
Tj et is the unit vector in the reconstructed jet direction and the exponent 0.7 was chosen 
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empirically. Then, following the method described in [44], the charge difference of the 
two di-jet pairs, 

AQ = (Qjeti + Qjet 2 ) — (Qjet s + Qjeu) , (8) 

was used to assign the charge of the individual W ± bosons. The more negative di-jet 
was tagged as originating from a W~ , and the other di-jet as the W + . The efficiency 
of this procedure was estimated from the simulation to be about 76% for events with 
correct jet pairing, using the minimal angle between the reconstructed di-jet system and 
the generated W boson to determine the correct pairing. As in the case of the jet pairing 
studies described above, the systematic uncertainty of this procedure was estimated by 
using the different parton shower and fragmentation models implemented in PYTHIA, 
HERWIG and ARIADNE. 

The ambiguity in the charge was taken into account by constructing a new variable: 

x q — P\v- (AQ) cos e w --(l- P w - (AQ)) cos 6 W - , (9) 

where cos 8w- is the polar angle of the di-jet pair assigned to the W~ and Pw- (AQ) is 
the probability that the di-jet pair originates from a W~. The value of Pw-(AQ) was 
obtained from the distribution of AQ in the simulated events. The couplings were then 
estimated from a binned extended maximum likelihood fit to the variable x q . 

5.3 Single W final state 

In the jjX final state, the couplings were extracted via a binned maximum likelihood 
fit to the distribution of the angle between the jets. This is a well-measured variable, and 
was found to be more sensitive to the coupling parameters than, say, the W production 
angle (a result which follows from the dynamics of the Feynman diagram (figure lc) pro- 
viding the dominant contribution to the sensitivity to the couplings in the jjX sample). 
The IX final state was analysed using a maximum likelihood fit to the number of events 
selected in the data, no further sub-division of the data being found to give a significant 
improvement to the experimental sensitivity. As mentioned in sections 4.4 and 4.5, the 
samples selected in these final states include contributions from some processes labelled 
as "background", but nonetheless with trilinear gauge couplings involved in their pro- 
duction mechanisms; in the fits performed, the relevant parameters were varied wherever 
they occurred in the production processes contributing to the events expected in the 
selected samples. 

Since only the WWj vertex occurs in the production of the Weu final state via a 
trilinear gauge coupling (as seen in figure 1(c)), the sensitivity of the single W channels 
to Agf is very poor, and fits to this parameter were not used in the results presented. The 
likelihood distributions from fits of the other two parameters, A 7 and Ak 7 , to the jjX 
and tX final states were combined, and the resulting distributions were subsequently 
combined with those from the jjtv and jjjj final states in the determination of the 
coupling parameters. 

6 Systematic uncertainties 

Sources of systematic uncertainty were considered which contribute to the results in 
all the final states analysed. Those arising in the analysis of the final states from WW 
production are described in section 6.1; the contribution to the total uncertainty from 
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each source to the results for each of the three coupling parameters determined from data 
in the jj£v and jjjj channels is given in tables 3 and 4, respectively. A similar study 
was performed for the couplings A 7 and Ak 7 determined from data in the single W final 
states. A summary is given in section 6.2 and the results are reported in table 5. 

6.1 WW final states 

W pair production cross-section and radiative corrections: 

The calculation of the W pair production cross-section was significantly improved in 
2000 [45] . The theoretical uncertainty in the relevant energy range was reduced from the 
level of 2% [46] quoted in previous publications [5] to 0.5% [19,47] via the inclusion of all 
0(a) corrections. The systematic uncertainties in the coupling parameters arising from 
this latest estimate of the uncertainty in the total WW cross-section are shown in the 
tables. 

The inclusion of the 0(a) corrections has also been shown to have a marked effect 
on the differential distributions [48], which could substantially affect the measurement 
of the gauge boson coupling parameters. The determination of the resulting systematic 
uncertainty in the determination of the couplings required the use of re-weighted events. 
The weights were generated using YFSWW in the simulation and were used according 
to the procedure described in [18]. The effect on the measurement of the couplings aris- 
ing from the theoretical uncertainty in the calculation of the radiative corrections was 
then obtained in two stages. First, one million fully simulated Monte Carlo events were 
produced at 189 GeV using the generators WPHACT, Racoon WW [47] and YFSWW. 
From a comparison of the couplings determined from analysis of these different samples, 
it was possible to estimate the systematic uncertainty from higher order electroweak cor- 
rections on the calculation of the Initial State Radiation. This was found to be negligible. 
Then, using the sample simulated with YFSWW, a comparison was made of two differ- 
ent Leading Pole approximation schemes, the so-called LPA-A and LPA-B schemes. The 
differences in the couplings determined from analysis of samples employing these two 
models was taken to represent the systematic error from the uncertainty of the depen- 
dence of the Double Pole Approximation on the assumed LPA scheme. It can be seen in 
the corresponding entries in tables 3 and 4 that this uncertainty gives rise to significant 
systematic errors in the measurement of the couplings. 

Background cross-sections and modelling: 

The theoretical uncertainty on the cross-sections of two- and four-fermion processes 
varies between 2% and 5%, depending on the process. A conservative estimate of the 
systematic error on the couplings was made by varying the predicted background cross- 
sections by ±5%. 

W mass and LEP beam energy: 

The systematic error arising from the uncertainty on the W mass used in the event 
simulation was evaluated using data samples generated with masses 1 GeV/c 2 above and 
below the nominal value. A linear interpolation was used to scale the systematic error to 
that which would arise from an uncertainty in the W mass of ±40 MeV/c 2 . 



13 



Source 




\ 

Ay 




WW cross-section 


0.0005 


0.0006 


0.007 


Radiative Corrections 


+0.005 
-0.002 


+0.004 
-0.002 


+0.012 
-0.015 


Background cross-section 


0.004 


0.003 


0.014 


W Mass 


0.001 


0.001 


0.002 


LEP beam energy 


0.0005 


0.0005 


0.001 


Luminosity 


0.0005 


0.0006 


0.007 


Fragmentation 


0.005 


0.005 


0.015 


Lepton tagging efficiency 


0.003 


0.003 


0.001 


Lepton charge assignment 0.005 


0.005 


0.003 


Jet reconstruction 


0.002 


0.002 


0.007 


Lepton reconstruction 


0.001 


0.001 


0.003 


Total 


+0.010 
-0.009 


+0.010 
-0.009 


+0.027 
-0.028 


Statistical errors 


+0.033 
-0.031 


+0.036 
-0.035 


+0.103 
-0.094 



Table 3: Contributions to the systematic errors on the couplings determined from data 
in the semi-leptonic final state, jjlv. Except where otherwise indicated, the errors are 
symmetric with respect to a change of sign of the parameters involved. The first 7 
sources listed in the table are considered to be fully correlated with the other channels. 
For comparison, the bottom row of the table lists the statistical errors on the couplings 
determined in the Optimal Variables analysis (also shown in table 6). 
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Source 








WW cross-section 


0.006 


0.008 


0.011 


Radiative corrections 


0.017 


0.016 


0.032 


Background cross-section 


0.003 


0.004 


0.009 


W mass 


0.003 


0.003 


0.005 


LEP beam energy 


0.001 


0.001 


0.001 


Luminosity 


0.006 


0.008 


0.011 


Fragmentation 


0.009 


0.012 


0.027 


Colour Reconnection 


0.008 


0.006 


0.012 


Bose Einstein 


0.002 


0.002 


0.005 


Simulation statistics 


0.008 


0.009 


0.012 


Selection efficiency 


0.005 


0.005 


0.007 


Event reconstruction 


0.004 


0.004 


0.008 


Total 


0.024 


0.025 


0.049 


Statistical errors 


+0.083 
-0.067 


+0.093 
-0.070 


+0.196 
-0.149 



Table 4: Contributions to the systematic errors on the couplings determined from data in 
the fully hadronic final state, jjjj. The first 7 sources listed in the table are considered 
to be fully correlated with the other channels. For comparison, the bottom row of the 
table lists the statistical errors on the couplings (also shown in table 6). 



Source 


A 7 


Ak 7 


Signal cross-section 


0.005 


0.037 


Background cross-section 


0.002 


0.002 


Selection efficiency 


0.011 


0.072 


Total 


0.011 


0.081 


Statistical errors 


+0.250 
-0.288 


+0.131 
-0.148 



Table 5: Contributions to the systematic errors on the couplings determined from data in 
the single W final states. For comparison, the bottom row of the table lists the statistical 
errors on the couplings (also shown in table 6). 
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The same method was used to estimate the systematic uncertainty due to the value 
of the Lep beam energy used in the simulation; samples were generated with differ- 
ent centre-of-mass energies and the errors were rescaled to the measured beam energy 
uncertainties [49]. 

Determination of the luminosity: 

The luminosity was determined from a measurement of Bhabha scattering and was 
affected by the experimental error on the acceptance (±0.5%) and the theoretical uncer- 
tainty on the cross-section (±0.12%) [50]. The estimated uncertainty on the luminosity 
was used to vary the normalisation of the simulation in the fits. 

Modelling of fragmentation: 

In order to assess the effect of the model used for the fragmentation of hadronic jets - 
JETSET final state QCD radiation and fragmentation, - correlated samples were analysed 
using the modelling of HERWIG and ARIADNE, and the differences in the fitted values 
of the coupling parameters noted. The largest discrepancies found were between JETSET 
and HERWIG and these were taken as a conservative estimate in each channel. 

Additional tests were performed in the fully hadronic final state using mixed Lorentz- 
boosted Z events [51], in which WW events are emulated using two events taken at the 
Z peak, and transforming them such that their superposition reflects that of a true fully 
hadronic WW event. These studies are also sensitive to systematic errors in the event 
reconstruction technique, and are discussed further in the relevant section below. 

Final state interactions: 

The measurement of the couplings in the fully hadronic final state is affected by final 
state interactions between the decay products of the two W bosons. Two effects were 
considered: the exchange of gluons between the quarks of different W bosons, known as 
Colour Reconnection, and Bose-Einstein correlations between pions. 

Colour Reconnection: 

In the reaction e + e~ — > W + W~ — > (qiq~2) (Q3Q4) the hadronisation models used in this 
analysis treat the colour singlets qiq 2 and q 3 q4 coming from each W boson independently. 
However, interconnection effects between the products of the two W bosons may be 
expected since the mean W lifetime is an order of magnitude smaller than the typical 
hadronisation times. This can lead to the exchange of coloured gluons between partons 
from the hadronic systems from different W bosons - the Colour Reconnection effect - 
in the development of the parton showers. This, in turn, can give rise to a distortion 
in the angular distributions of the final hadronic systems used to estimate the primary 
quark directions in the determination of the triple gauge coupling parameters from jjjj 
data. These effects can be large at hadronisation level, due to the large numbers of 
soft gluons sharing the space-time region, and have been studied by introducing colour 
reconnection effects into various hadronisation models. The most studied model is the 
Sjostrand-Khoze "Type 1" model (SKI) [52], and this was used for the evaluation of the 
systematic uncertainty in the analysis reported here. The model is based on the Lund 
string fragmentation phenomenology, in which the volume of overlap between two strings, 
and hence the colour reconnection probability, is represented by a parameter, k. 
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In this paper, the systematic uncertainty was estimated using the SKI model with 
full colour reconnection (i.e. k = oo). This is a highly conservative assumption when 
compared with the direct measurements of colour reconnection reported by DELPHI [53] 
and by other LEP experiments [54-56]. Symmetric systematic errors were applied to 
the gauge coupling parameters, representing the difference observed between full colour 
reconnection and no effect from this source. 

Bose-Einstein correlations: 

Correlations between final state hadronic particles are dominated by Bose-Einstein 
correlations, a quantum mechanical effect which enhances the production of identical 
bosons close in phase space. The net effect is that multiplets of identical bosons are pro- 
duced with smaller energy-momentum differences than non-identical ones. This, again, 
can affect the estimation of the primary quark directions in data from hadronically de- 
caying W bosons. Bose-Einstein correlations between particles produced from the same 
W boson affect the normal fragmentation and are treated implicitly in the fragmenta- 
tion uncertainties constrained by data from Z decays, while correlations between pairs 
of particles coming from different W bosons cannot be constrained or safely predicted 
by the information from single hadronically decaying vector bosons, and are estimated in 
various models. We have used the LUBOEI BE 32 model [57] to estimate the systematic 
uncertainty in the determination of gauge coupling parameters from the present data. In 
this model, Bose-Einstein correlations are described using two parameters: the correla- 
tion strength, A, and the effective source radius, R. Applying the model with parameters 
A = 1.35 and R = 0.6 fm, symmetric systematic errors on the gauge coupling parame- 
ters were estimated by taking the difference between the values obtained assuming the 
presence of Bose-Einstein correlations only within each W and those obtained assuming 
correlations both within and between W bosons. Taking into account the reported results 
of measurements of Bose-Einstein correlations by Delphi [58] and in other LEP [59-61] 
experiments, this again represents a conservative estimate of the effect from this source. 

Statistics of simulated samples and selection efficiency: 

The statistical error on the number of simulated events assigned to each data bin 
was convoluted in the fitting method for fits to the data in the semi-leptonic channel; 
the fitting method ensures that this systematic error is negligible with the large statis- 
tics available. In the fully hadronic channel, the distribution of simulated events used 
in the binned extended maximum likelihood fit was varied according to the statistical 
uncertainties of the bin contents. 

The uncertainty due to the event selection efficiency was used to vary the normalisation 
of the simulation in the fits. 

Lepton tagging efficiency and charge assignment: 

Comparisons were made between fully simulated events and real Z events to estimate 
the possibility of having different lepton tagging efficiencies in the data. The systematic 
uncertainty was estimated assuming 1% mis-tagging for muons and for electrons in the 
barrel region and 5% for electrons in the forward region of the detector. The value shown 
in table 3 represents the combined effect from both lepton types, with the dominant 
contribution coming from mis-tagged electrons. However, the effect is reduced as mis- 
tagged electrons or muons can be retrieved by the single prong tau selection. 
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The effect of wrongly assigned lepton charge was estimated using data simulated at the 
Z pole by counting the numbers of di-lepton events in which the two leptons are assigned 
the same charge. A mis-assignment rate of 0.1% was found for all lepton candidates 
except for electrons in the forward region, where the rate rose to 6%. The systematic 
error was calculated by randomly changing the charge of the lepton candidate in the fits 
with these probabilities, and the value shown in table 3 shows the combined effect of 
these assumed uncertainties. 

Event reconstruction: 

The effect of possible systematic errors in the event reconstruction technique was 
estimated using comparisons between data and simulation. This was performed in two 
ways: firstly, by comparing significant variables used in the analysis in data and simulation 
and computing the effect of the discrepancy seen; and secondly, by directly computing 
changes in the results using mixed Lorentz-boosted Z (MLBZ) events, mentioned above 
in the section on systematic errors resulting from the modelling of fragmentation. 

In the semi-leptonic channel, the systematic uncertainty in the couplings due to un- 
certainties in the lepton and jet energies and angular distributions was estimated using 
comparisons between data and simulated events at the Z peak. The estimated uncertain- 
ties on the jet energies and angles were found to be 5% and 7.5 mrad, respectively. The 
uncertainty on the muon momentum was found to be 1%, while for electron momenta 
uncertainties of 1% and 5% were estimated in the barrel and forward regions, respectively. 
Appropriate smearings were applied to these resolutions in the simulation of jjiu events 
and the resulting shifts in the values of the couplings were taken to be the systematic 
uncertainties. They are reported in table 3 as the systematic errors arising from jet and 
lepton reconstruction. 

In the fully-hadronic channel, the uncertainties in the event reconstruction were es- 
timated using MLBZ events from both real and simulated data at the Z peak. As 
described above in the discussion of the modelling of quark fragmentation, the MLBZ 
method emulates WW events using two events taken at the Z peak, rotating them and 
Lorentz-boosting them such that their superposition reflects that of a true WW event. 
The detector effects are thus included in as realistic a manner as possible. In order to 
estimate these effects on the determination of gauge coupling parameters in WW^t jjjj 
events, the ratio of selection efficiencies, r, of MLBZ data events to MLBZ simulated 
events was determined as a function of the simulated W production angle. The ratio 
r (cos 6w-) was then applied to simulated WW samples and the gauge coupling anal- 
ysis described in section 5.2, which uses the W production angle, was repeated. The 
differences between the results with and without application of the ratio were taken as 
systematic errors and are reported in table 4. The systematic uncertainty evaluated 
by this method represents a conservative estimate, as it includes both the inaccuracies 
in the modelling of detector effects and most of the deviations induced by the applied 
fragmentation model. 

An additional problem, not included in the effects considered above, has been en- 
countered in the reconstruction of charged tracks in the forward region of Delphi [51], 
leading to a small error in the reconstructed direction of forward tracks in both simulated 
and real data. Its effects were shown to be negligible in a previous Delphi analysis [6] 
involving fits to binned data of production and decay distributions in WW production, 
and, in a study of the current data in the jjlv final state at 200 GeV, have also been 
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found to be negligible in comparison to the other correlated systematic errors considered. 
No systematic error has therefore been included from this source. 

6.2 Single W final states 

Systematic errors arising from the uncertainty in the signal (Wei/) cross-section were 
estimated by varying the cross-section by ±5% and noting the effect on the fitted coupling 
parameters. Similarly, cross-sections of other contributing channels were varied by ±2%, 
and the fits repeated. The maximum changes to the fitted parameters in the combined 
jjX and £X final states were taken as systematic errors, and are reported in table 5 
as the contributions from signal and background cross-sections, respectively. Systematic 
errors arising from the uncertainty in the selection efficiency were estimated from the 
statistical errors in the numbers of simulated events, and are also reported in the table. 
No other significant source of systematic error was identified in these channels. 

7 Results 

The procedure used to combine the results from the three channels and the results 
obtained are described in the following sections. 

7.1 Combination procedure 

The combination was based on the individual likelihood functions from the samples 
of the three final states, jj£v, jjjj and Weu, included in the analysis. Each final state 
sample provides the negative log likelihood, -ln£, at each centre-of-mass energy, as a 
function of the coupling parameters for inclusion in the combination. 

The ln£ functions from each channel include statistical errors as well as those sys- 
tematic uncertainties which are considered as uncorrelated between channels. For both 
single- and multi-parameter combinations, the individual ln£ functions were combined. 
It is necessary to use the In £ functions directly in the combination, since in some cases 
they are not parabolic, as discussed extensively in [62], and hence it is not possible to 
combine the results properly by simply taking weighted averages of the measurements. 

The following sources of systematic uncertainty were assumed to be correlated between 
the semi-leptonic and fully hadronic channels: WW cross-section, radiative corrections, 
background cross-section, W mass, beam energy, luminosity and fragmentation. The pro- 
cedure used was based on the introduction of an additional free parameter to take into 
account each correlated source of systematic uncertainty. These parameters are treated 
as shifts on the fitted parameter values, and are assumed to have Gaussian distribu- 
tions. A simultaneous minimisation of both sets of parameters (coupling parameters and 
systematic uncertainties) was performed on the log-likelihood function. 

In detail, the combination proceeded in the following way: the set of measurements 
from the three channels jjlv, jjjj and single W is given with statistical plus uncor- 
related systematic uncertainties in terms of likelihood curves — \aCl^(x), — h\C q ™(x) 
and — In C s s l t ^f eW (x) , respectively, where x is the coupling parameter in question. Also 
given are the shifts for each of the totally correlated sources of uncertainty mentioned 
above, each source S giving rise to systematic errors a^ lv and cr^ ■ Additional parame- 
ters A s are then included in the likelihood sum in order to take into account a Gaussian 
distribution for each of the systematic uncertainties. The procedure then consisted in 
minimising the function 
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- In C total = - £ In C c stat {x - J2(^ S )) + E V~ ' (10) 

C 5 5 

where x and the A 5 are the free parameters, the sum over C runs over the three channels 
and the sum over S runs over the seven sources of correlated systematic uncertainty. 
The resulting uncertainty on x takes into account all sources of uncertainty, yielding 
a measurement of the coupling with a precision which includes the errors from both 
statistical and systematic sources. The projection of the minima of the log-likelihood as 
a function of x gives the combined log-likelihood curve including statistical and systematic 
uncertainties. 



7.2 Results 

The data taken by Delphi between 1998 and 2000 were collected at centre-of-mass 
energies between 189 and 209 GeV. The results for the measurement of the couplings 
from single parameter fits to the data in the different channels are given in table 6 and 
the likelihood curves from these fits are shown in figure 9. The results of the simultaneous 
fits to the data for all combinations of two parameters (Agf-A 7 , Agf-An^ and A 7 -Ak 7 ) 
are given in table 7. The corresponding likelihood contours are shown in figure 10. The 
result from the simultaneous fit to all three couplings is given in table 8 and the likelihood 
contours corresponding to the intersections of the three 2-parameter planes containing 
the minimum of the distribution with the three-dimensional 3-parameter likelihood dis- 
tribution are shown in figure 11. 

It may be noted from the results shown in tables 6-8 that the 68% and 95% confidence 
levels obtained in the 3-parameter fit are somewhat narrower than those obtained in the 
1-parameter fit to the same parameter. This is not expected if the likelihood distributions 
are strictly Gaussian in form. However, such an effect is also observed in analysis of a 
significant fraction (5%) of simulated event samples of the same size as the experimental 
sample. As has been pointed out in previous studies of both simulated [63] and experi- 
mental [10] samples, the quadratic dependence of the cross-section on the couplings we 
consider does indeed lead to non-Gaussian likelihood distributions, which can thus ex- 
plain this behaviour. The results for the multidimensional fits, in particular those for the 
3-parameter fit, should therefore be viewed with this constraint on their interpretation 
in mind. 

The result from the simultaneous fit to A 7 and Ak 7 can be converted to a measurement 
of the magnetic dipole moment, /j,w, and the electric quadrupole moment, qw, of the W + 
boson using the relations given in equations (2) and (3). The resulting two-parameter fit 
is shown in figure 12. The fitted values of /iw and qw are 

Hw ■ 2m w /e = 2.027+°;°^ 
and q w -m 2 w /e = -1.025^ , 

where the errors include both statistical and systematic contributions. These results may 
be compared with the Standard Model predictions of 2 and —1 for these two quantities, 
respectively. 
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Channel Agf A 7 Ak 



7 



jj£v (Optimal Variables) -0.0241°;°^ 0.006^;°^ 0.014+°;^ 
jjiu (Angular Variables) 0.0061^39 0.019l°;°43 -0.091+ 096 



0.085 

nn -o.o30+Hf 7 -o.032lg:8?g o.osil^ 

single W 0.037^1° 0.0271°;^ 

Combined -0.025+°;°^ 0.0021°;°;* 0.024l°;°£ 



Table 6: The results for single parameter fits to the couplings in the individual channels. 
In each fit, the other two couplings were held at their Standard Model values. The errors 
given for the individual analyses are statistical; the systematic contributions are given in 
tables 3, 4 and 5. As indicated in the text, the Angular Variables analysis of the jjiv 
final state was performed as a cross-check, the values in the combination of all three 
channels being obtained using the results from the jjlv Optimal Variables analysis. The 
combined results also contain the systematic errors, included via the combination method 
described in the text. 



Parameter 


68% C.L. 


95% C.L. 


Correlations 


Agf 

A 7 


U.U4O_ 040 

n nQ7+o.o45 

U.UO 1 _ .044 


[-0.123, +0.035] 
[-0.051, +0.124] 


1.0 -0.49 
-0.49 1.0 


Aj/f 
Akj 


-0.033^1 
0.0591°;°?* 


[-0.097, + 0.032] 
[-0.093, + 0.233] 


1.0 -0.41 
-0.41 1.0 


A-y 

Ak 7 


-0.002^1 
0.0281™?? 


[-0.070, + 0.067] 
[-0.120, +0.198] 


1.0 0.10 
0.10 1.0 



Table 7: The measured central values, one standard deviation errors and limits at 95% 
confidence level obtained by combining the different channels in the 3 two-parameter 
fits. Since the shape of the log-likelihood is not parabolic, there is some ambiguity in 
the definition of the correlation coefficients and the values quoted here are approximate. 
In each fit, the listed parameters were varied while the remaining one was fixed to its 
Standard Model value. Both statistical and systematic errors are included. Note that the 
68% and 95% confidence limits reported here refer to single-parameter errors (in contrast 
to those shown in the two-parameter plots of figure 10) and are defined by A In C = +0.5 
and Aln£ = +1.92, respectively. 
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Parameter 


68% C.L. 


95% C.L. 


Correlations 
Agf A 7 A/t 7 


Ak 7 


-o-o6o±8:8|S 

0.038lffi 

o.o77±g : g?s 


[-0.118, +0.002] 
[-0.027, + 0.099] 
[-0.050, +0.218] 


1.0 -0.55 -0.41 
-0.55 1.0 -0.04 
-0.41 -0.04 1.0 



Table 8: The measured central values, one standard deviation errors and limits at 95% 
confidence level, obtained by combining the different channels in the three-parameter 
fit. Since the shape of the log-likelihood is not parabolic, there is some ambiguity in the 
definition of the correlation coefficients and the values quoted here are approximate. Both 
statistical and systematic errors are included. Note that the 68% and 95% confidence 
limits reported refer to single-parameter errors and are defined by Aln£ = +0.5 and 
Aln£ = +1.92, respectively. 

8 Conclusions 

The data taken by Delphi at centre-of-mass energies between 189 and 209 GeV have 
been used to probe the non-Abelian nature of the Standard Model. Limits have been 
placed on the trilinear gauge boson couplings which describe the WW Z and WW'j ver- 
tices; in particular, reactions leading to W pair production and single W production have 
been used to set limits on the parameters Agf , A 7 and A/t 7 . The combined results for 
fits to a single parameter, where the other two parameters were held at their Standard 
Model values, are: 

-0.084 < Agf < 0.039 , 
-0.065 < A 7 < 0.071 , 
and -0.129 < Ak 7 < 0.182 

at 95% confidence level. Fits were also made where two or three parameters were allowed 
to vary simultaneously. No deviations from the Standard Model predictions have been 
observed. 
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Figure 2: Distribution of the Optimal Variables oo\ and c^ 1 (the coefficients, respectively, 
of Agf and of {Agf ) 2 in the expansion of the differential cross-section in terms of Optimal 
Variables) for semi-leptonic data at 200 GeV. The points represent the real data, the solid 
lines the expected distributions for the SM value of the coupling, and the dashed lines 
the expected distributions for the non-SM values Agf=±l. The shaded area represents 
the background. The simulated distributions are normalised to the same luminosity as 
the data. 
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Figure 3: Distribution of the Optimal Variables uj and w 22 (the coefficients, respectively, 
of A 7 and of (A 7 ) 2 in the expansion of the differential cross-section in terms of Optimal 
Variables) for semi-leptonic data at 200 GeV. The points represent the real data, the solid 
lines the expected distributions for the SM value of the coupling, and the dashed lines 
the expected distributions for the non-SM values A 7 =±l. The shaded area represents 
the background. The simulated distributions are normalised to the same luminosity as 
the data. 
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Figure 4: Distribution of the Optimal Variables oo\ and cuf 3 (the coefficients, respectively, 
of Ak 7 and of (Ak 7 ) 2 in the expansion of the differential cross-section in terms of Optimal 
Variables) for semi-leptonic data at 200 GeV. The points represent the real data, the solid 
lines the expected distributions for the SM value of the coupling, and the dashed lines 
the expected distributions for the non-SM values A/t 7 = ±3. The shaded area represents 
the background. The simulated distributions are normalised to the same luminosity as 
the data. 
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Figure 5: Distribution of the Optimal Variables u^ 2 , u 23 and u^ 3 (the coefficients, respec- 
tively, of Agf -A 7 , A 7 -A/t 7 and Agf -Ak 7 in the expansion of the differential cross-section 
in terms of Optimal Variables) for semi-leptonic data at 200 GeV. The points represent 
the real data, the solid lines the expected distributions for SM values of the couplings, and 
the dashed lines the expected distributions for the non-SM values of the couplings shown 
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Figure 6: The distribution of cos Oyy+, the cosine of the polar angle of the W + in semi-lep- 
tonic events, at a centre-of-mass energy of 200 GeV. All three plots show the data (rep- 
resented by points), the Standard Model prediction (the solid line) and the predicted 
background (the darker shaded region). Each plot also shows predictions for non-Stan- 
dard Model values of a coupling a: in a), a = Agf , in b) a = A 7 , and in c) a = Ak 7 . 
The simulated distributions are normalised to the same luminosity as the data. 
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Figure 7: The distribution of cos 6i, the cosine of the polar angle of the lepton in semi-lep- 
tonic events with respect to the incoming e ± of the opposite sign, at a centre-of-mass 
energy of 200 GeV. All three plots show the data (represented by points), the Standard 
Model prediction (the solid line) and the predicted background (the darker shaded re- 
gion). Each plot also shows predictions for non-Standard Model values of a coupling a: 
in a), a = Agf, in b) a = A 7 , and in c) a = Ak 7 . The simulated distributions are 
normalised to the same luminosity as the data. 
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Figure 8: The distribution of cos Bxw, the cosine of the angle between the directions of the 
lepton and the hadronic W in semi-leptonic events, at a centre-of-mass energy of 200 GeV. 
All three plots show the data (represented by points), the Standard Model prediction (the 
solid line) and the predicted background (the darker shaded region). Each plot also shows 
predictions for non-Standard Model values of a coupling a: in a), a = Agf , in b) a = A 7 , 
and in c) a = Ak 7 . The simulated distributions are normalised to the same luminosity 
as the data. 
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Figure 9: The log-likelihood curves from single parameter fits to the data, combining 
results from the semi-leptonic, fully hadronic and single W final states. The curves include 
contributions from both statistical and systematic effects 
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Figure 10: The log-likelihood contours for two-parameter fits to the data, combining 
results from the semi-leptonic, fully hadronic and single W final states. The plots include 
contributions from both statistical and systematic effects 
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Figure 11: Intersections of the 68% and 95% confidence level 3-parameter log-likelihood 
surfaces with the three 2-parameter planes containing the minumum of the 3-parameter 
likelihood fit. The figures combine the results from the semi-leptonic, fully hadronic and 
single W final states and include contributions from both statistical and systematic effects 
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Figure 12: The log-likelihood contours for a two-parameter fit to qw and fiw, respectively 
the electric quadrupole and magnetic dipole moments of the W + boson, obtained from the 
simultaneous fit to A 7 and Ak t Results from the semi-leptonic, fully hadronic and single 
W final states have been combined in the plot, and contributions from both statistical 
and systematic effects are included. The Standard Model expectations for the quantities 
plotted are: qwrn^/e = —1 and fiw^mw/e = 2. 



